Odd-frequency pairing and Josephson effect at superconducting interfaces 
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We demonstrate that, quite generally, the spin-singlet even-parity (spin-triplet odd-parity) pair potential in a 
superconductor induces the odd-frequency pairing component with spin-singlet odd-parity (spin-triplet even- 
parity) near interfaces. The magnitude of the induced odd-frequency component is enhanced in the presence 
of the midgap Andreev resonant state due to the sign change of the anisotropic pair potential at the interface. 
The Josephson effect should therefore occur between odd- and even-frequency superconductors, contrary to the 
standard wisdom. A method to probe the odd-frequency superconductors is proposed. 

PACS numbers: 74.45,+c, 74.50.+r, 74.20.Rp 



It is well known that Josephson coupling occurs between 
superconductors with the same order parameter symmetry fl]] . 
Generally, symmetries with respect to momentum, spin and 
time are considered. On the basis of symmetry with respect to 
time, two classes of superconductors are introduced, referred 
to as odd-frequency and even-frequency superconductors. In 
accordance with the Pauli principle, the even-frequency state 
is characterized by spin-singlet even-parity or spin-triplet odd- 
parity order parameter, while odd-frequency superconductors 
belong to spin-singlet odd-parity or spin-triplet even-parity 
pairing state. It was suggested in Ref. yfl from basic symme- 
try arguments that the first order Josephson coupling should 
be absent between odd- and even-frequency superconductors. 
However, as shown in the present paper, an odd-frequency 
component of the pair potential is quite generally induced 
at interfaces in superconductors with even-frequency bulk 
pair potential, while an even-frequency component is induced 
at interfaces in odd-frequency superconductors. Therefore 
Josephson coupling between even- and odd-frequency super- 
conductors should be possible. 

Up to now, almost all of known superconductors belong to 
the symmetry class of the even-frequency pairing. The possi- 
bility of the odd-frequency pairing state in a uniform system 
was proposed for 3 He in y|] and for a superconducting state 
involving strong correlations J2,|3,|5l0]. Recently, the realiza- 
tion of the odd-frequency pairing state in inhomogeneous sys- 
tems was proposed by Ref. 7 in ferromagnet/superconductor 
heterostructures. This issue was further addressed in several 
related studies |8|,|9|]. Recent experiments provided evidence 
for such anomalous pairing states [10]. Furthermore, two of 
the present authors predicted that the odd-frequency pairing 
state can be induced in a diffusive normal metal attached to a 
spin-triplet superconductor . However, in these examples 
the odd-frequency pairing is realized due to spin-triplet cor- 



relations. A question naturally arises whether the spin-triplet 
ordering is a necessary condition for the existence of the odd- 
frequency pairing state at superconducting interfaces. In order 
to address this issue, one should consider self-consistently the 
pair potential near the interface. 

The classic example of an inhomogeneous superconduct- 
ing system is a ballistic normal metal/superconductor (N/S) 
junction with an even-frequency superconductor. Particularly 
interesting are the cases when S is a spin-singlet <i-wave or a 
spin-triplet p-wave superconductor since the sign change of 
the pair potential probed by quasiparticles injected and re- 
flected by the N/S interface can occur in such junctions. As 
a result, the pair potential is suppressed near the interface 
1 1211 and the midgap Andreev resonant states (MARS) are 
formed [13]. The appearance of unusual charge transport in 
the presence of MARS 111 311 suggests the presence of underly- 
ing anomalous pairing states. 

In this Letter we address the issue of odd-frequency pairing 
in the generic case of a N/S interface, where superconductor 
S has an even-frequency pairing state in the bulk. We will 
use the quasiclassical Green's function formalism where the 
spatial dependence of the pair potential is determined self- 
consistently. We will show that, quite generally, the odd- 
frequency component is induced near superconducting inter- 
faces due to the spatial variation of the pair potential [i4J. If 
a superconductor has an ESE or ETO pairing state in the bulk, 
the order parameter at the interface has respectively an odd- 
frequency spin-singlet odd-parity (OSO) or an odd-frequency 
spin-triplet even-parity (OTE). In the absence of MARS (e.g. 
ESE superconductor with s-wave or d x 2 ^ y 2 -wave symmetry 
of the order parameter), the magnitude of the odd-frequency 
component of the pair amplitude decreases with a decrease of 
the transmission coefficient T m across the interface. On the 
other hand, in the presence of MARS (e.g. ETO superconduc- 
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tor with p^-wave or ESE one with d^-wave symmetry of the 
order parameter), the magnitude of the odd-frequency compo- 
nent is enhanced with a decrease of T m . Similarly, in the case 
of a bulk odd-frequency superconductor, an even-frequency 
component of the order parameter is induced at the interface. 
An important application of the above results is the existence 
of Josephson coupling between bulk odd- and even-frequency 
superconductors to the first order in T m . 

In the following, we consider an N/S junction as the sim- 
plest example of a non-uniform superconducting system with- 
out impurity scattering. Both the ESE and ETO states are 
considered in the superconductor. As regards the spin-triplet 
superconductor, we choose S z = for the simplicity. We 
assume a thin insulating barrier located at the N/S interface 
(x = 0) with N (x < 0) and S(x > 0) modeled by a delta 
function HS(x), where H is the magnitude of the strength of 
the delta function potential. The reflection coefficient of the 
junction for the quasiparticle for the injection angle is given 
by R = Z 2 /{Z 2 + 4 cos 2 0) = 1 - T m with Z = 2H/v F , 
where (— tt/2 < 9 < tt/2) is measured from the normal to 
the interface and vf is the Fermi velocity. The quasiclassi- 
cal Green's function in a superconductor is parameterized in 
terms of 

9± = fi±h + /2±t 2 + g±r 3 , g± = 1 (1) 

with Pauli matrices fj (i = 1 — 3) and unit matrix 1. Here, 
the index + (— ) denotes the left (right) going quasiparticles 
lfl6ll . It is possible to express the above Green's function as 
/i± = ±i(F±+D±)/(l-D ± F ± ), h± = -(F±-D±)/(l- 
D±F±),wdg± = (1 + D±F±)/(1-D±F±), where D± and 
F± satisfy the Eilenberger equations in the Riccatti parame- 
terization Il7ll 

v Fx d x D± = -A±(x)(l-D 2 ± ) + 2uj n D± (2) 

v Fx d x F ± = -A±(x)(l - Fl) - 2uo n F±. (3) 

where vf x is the x component of the Fermi velocity, uj n = 
2ixT[n + 1/2) is the Matsubara frequency, with tempera- 
ture T. A+(x) (A_(x)) is the effective pair potential for 
left going (right going) quasiparticles. Since the interface 
is flat, F± = -RD T holds at a; = lfl7h. Here we con- 
sider the situation without mixing of different symmetry chan- 
nels for the pair potential. Then A±(x) is expressed by 
A±(x) = A(x)$±(0)Q(x) with the form factor $±(9) given 
by 1, cos 29, ± sin 29, ±cos0, and sin0 for s, d x 2_ y 2 , d xy , 
p x andp a -wave superconductors, respectively. A(x) is deter- 
mined by the self-consistency equation 

ot r 7r / 2 

A(x) = V / d9G{9)h+ (4) 

n>l * 2 

with G(9) = 1 for the s-wave case and G(0) = 2<5>(6) for 
other cases, respectively 111 811 . Tc is the transition temperature 



of the superconductor. The condition in the bulk is A (00) = 
A. Since the pair potential A(x) is a real number, the resulting 
fi± is an imaginary and / 2 ± is a real number. 

In the following, we explicitly write fi± = fi±(oj n ,9), 
h± = /a±K,0), F± = F±{iu n ,9) and D ± = D± (w„,6>). 
For x — 00, we obtain fx±(u> n ,6) — and f2±(u> n ,8) = 
A$ ± (9) Note that fi+(uj n , 9) becomes finite due to 

the spatial variation of the pair potential and does not exist as 
the bulk. From eqs. (f2]) and (01, we can show D±(— u> n , 9) = 
l/D ± (oj n ,6) and F±(-u n ,0) = 1/F±(u n ,6). After a sim- 
ple manipulation, we obtain fi±(u> n ,9) = — /i±(— 
and f2±(u n ,9) = /a±(— w n ,^) for any x. It is remarkable 
that functions fi±(ui n ,0) and fa±(u> n , 9) correspond to an 
odd-frequency and an even-frequency one of the pair ampli- 
tude, respectively fl. 

Next, we shall discuss the parity of these pair ampli- 
tudes. The even-parity (odd-parity) pair amplitude should sat- 
isfy the relation fi±(u n ,6) = f iT (u n ,-0) (fi±(ui n ,9) = 
—fi^(u) n ,—8)), with i = 1,2. For an even(odd)-parity su- 
perconductor $±(-0) = $ T (0) ($±(-0) = -$ T (0)). 
Then, we can show that D±(-0) = D T {6) and F±(-0) = 
F T (9) for even-parity case and D±{—6) = —D T (9) and 
F±{— 0) = — F T (9) for odd-parity case, respectively. The 
resulting fi±(u) n ,6) and / 2 ±(w„,0) satisfy fi±{u n ,9) = 
-fiT&n,-®) and / 2 ±(w rl ,0) = f 2T (u n ,~9) for an even- 
parity superconductor and fi±(u> n ,6) — f\^(uj n , — 0) and 
/2±(w„, 0) = —f2zfz(w n > —9) f° r an odd-parity superconduc- 
tor. Note that the parity of the odd-frequency component 
fi±(u n , 9) is different from that of the bulk superconductor 
for all cases. 

Let us now focus on the values of the pair amplitudes at 
the interface x — 0. We concentrate on two extreme cases 
with (1) $+(0) = $_(0) and (2) $+(0) = -$-(0). In 
the first case, MARS is absent since there is no sign change 
of the pair potential felt by the quasiparticle at the interface. 
Then D + = D- is satisfied. On the other hand, in the second 
case, MARS is generated near the interface due to the sign 
change of the pair potential. Then D + = —D- is satisfied 
JI3I1 . At the interface, it is easy to show that fi± = ±i(l — 
R)D+/{1 + RD\) and f 2± = (1 + R)D+/(1 + RD 2 + ) for 
Case (1) and f 1± = i{\ + R)D+/(1 - RD\) and f 2 ± = 
±(1 - R)D+/{1 - RD 2 + ) for Case (2), respectively, where 
the real number D+ satisfies | D + | < 1 for u n 7^ 0. For Case 
(1), the magnitude of fi± is always smaller than that of / 2 ±. 
For Case (2), the situation is reversed. In the low transparent 
limit with R — > 0, only the f\± is nonzero. Namely, only the 
even-frequency (odd-frequency) pair amplitude exists at the 
interface without (with) sign change of the pair potential. 

In order to understand the angular dependence of the pair 
amplitude in a more detail, we define fx and / 2 for — tt/2 < 
9 < 3vr/2 with / 1(2) = A (2)+ (0) for -vr/2 < < vr/2 and 
A(2) — /i( 2 )-( 7r — 0) f° r 7r /2 < < 37r/2. We decompose 
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/i(2) into various angular momentum component as follows, 



'£ (a)s-wave junction 



1(2) 



E 



S^ 2)) sm[m6} + ^ C^ 1(2)) cos[to0] (5) 



with m = 21 + 1 for odd-parity case and m = 21 for even 
parity case with integer I > 0, where I is the quantum number 
of the angular momentum. It is straightforward to show that 



the only nonzero components are (1) C, 



(2) 
21 



or cL. 



,,2-wave. 



and cff for p^-wave, and (4) S^m an< ^ ^21' f° r £V wave 
junctions, respectively. 

Below we illustrate the above results by numerical calcula- 
tions. As typical examples, we choose s-wave and p^-wave 
pair potentials. Although both /1 and f% have many compo- 
nents with different angular momenta, we focus on the lowest 



(2) Sg> 



and S. 



(i) 



and C^ +l for s 
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(iu} n ,x) = C ( 
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C^, O s (iuj n ,x) = Cq 1 ', and O px (iLu n ,x) — and 
choose ?'oj„ = iixT with temperature T = 0.05Tc. For the 
,s-wave case, the pair potential is suppressed only for high 
transparent junctions (see Fig. la). The odd-frequency com- 
ponent O px (iirT,x) is enhanced for Z = near the inter- 
face where the pair potential A(x) is suppressed, while the 
even-frequency component E s (iirT,x) remains almost con- 
stant in this case. For low transparent junctions, the magnitude 
of O px (iirT,x) is negligible (see Fig. lb). For the p^-wave 
junction with Z = 0, although the odd-frequency component 
O s (iirT, x) is enhanced near the interface, it is smaller than 
the even-frequency one E px (iirT, x) (see Fig. lc). For the low 
transparent junction, the magnitude of O s (iirT, x) is strongly 
enhanced near the interface and becomes much larger than 
the magnitude of E px (iirT, x) (see Fig. Id). It is a remarkable 
fact that the MARS is reinterpreted as the manifestation of the 
odd-frequency pair amplitude at the interface. 

If we choose d^-wave pair potential, the magnitude of 
the odd-frequency component is enhanced at the interface 
for large Z. In this case, the OSO state is generated while 
the OTE state is generated for p^-wave case. The s-wave 
isotropic component that is robust against the impurity scat- 
tering Ql UM appears only in the p^-wave case. The OSO 
state cannot penetrate into a diffusive normal metal while the 
OTE state can. Thus we can naturally understand the presence 
of proximity effect with MARS in p^-wave junctions and its 
absence in dry-wave junctions 111 1 , 



In the following, we discuss the manifestations of the odd- 
frequency pairing in the Josephson effect between even- and 
odd-frequency superconductors. For this purpose, we first ex- 
tend the above discussion regarding Eqs. (l)-(5) to the case 
of bulk odd-frequency pairing in a superconductor. In this 
case, A±(x) depends on a Matsubara frequency A±(x) — 
A±(x, Lo n ) with A±(x, oj n ) — — A±(x, — io n ). Onecanshow 
that /2± (/i±) becomes an odd-frequency (even-frequency) 
pair amplitude and ETO (ESE) pairing state is induced at the 
interface for the OTE (OSO) bulk superconductor. For the low 
transparency limit, only odd-frequency states exist at the inter- 
face for A±(x,oj n ) = A T (x,uj n ) where the sign change of 




(b)s-wave junction 
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FIG. 1: Spatial dependence of the normalized pair potential (solid 
line) even-frequency pair amplitude (dotted line) and odd-frequency 
pair amplitude (dashed line). For a bulk state, we choose ESE s- 
wave and ETO p^-wave superconductor, (a) and (c): fully transpar- 
ent junctions with Z = 0. (b) and (d): low transparent junctions 
with Z = 5. E s (iirT, x) and E Vx (iirT, x) are the even-frequency 
components of the pair amplitude, O s (iirT, x) and O px (*7rT, x) are 
the odd-frequency components. The distance x is normalized by the 
v F /A. 



the pair potential is absent. On the other hand, in the pres- 
ence of the sign change, i.e., A±(x,u n ) — —A T (x,u> n ), 
only even-frequency states exist. 

To summarize the above results, we present the relation be- 
tween the order parameter symmetry in the bulk superconduc- 
tor and at the interface for T m — * (free surface). As shown 
in the Table I below, there are eight distinct cases which corre- 
spond to different combinations of the bulk pairing symmetry 
and the behavior of the orbital part of the bulk pair potential 
with respect to reflection from the interface. 



Let us discuss the Josephson coupling at the interface be- 
tween even-frequency and odd-frequency superconductors to 
the first order in the interface transparency coefficient T m , as- 
suming that spin-flip scattering at the interface is absent. Ac- 
cording to the Table I, these are 16 possible combinations of 
pairing symmetries in two superconductors. Due to the dif- 
ference of the spin structure of Cooper pairs, the Josephson 
coupling is absent for the following combinations: (l)-(7), 
(l)-(8), (2)-(7), (2)-(8), (3)-(5), (3)-(6), (4)-(5), and (4)-(6). 
The Josephson coupling is also absent for the combinations 
(l)-(5), (2)-(6), (3)-(7) and (4)-(8), since the odd- and even- 
frequency pairing states are realized on both sides of the in- 
terface. This result is consistent with the previous prediction 
021]. The remaining four combinations (l)-(6), (2)-(5), (3)-(8), 
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bulk state 


sign change 


interface state 


(1) 


ESE (s or d x 2 _ y i -wave) 


No 


ESE 


(2) 


ESE (d IH -wave) 


Yes 


OSO 


(3) 


ETO (py-wave) 


No 


ETO 


(4) 


ETO (p^-wave) 


Yes 


OTE 


(5) 


OSO (p^-wave) 


No 


OSO 


(6) 


OSO (pn-wave) 


Yes 


ESE 


(7) 


OTE (s or d x 2 _ v i -wave) 


No 


OTE 


(8) 


OTE (d IH -wave) 


Yes 


ETO 



TABLE I: The relation between the symmetry of the bulk supercon- 
ductor and that of the pair amplitude at the interface in the low trans- 
parent limit. The allowed symmetry of the Cooper pair in accordance 
with Pauli's rule are even-frequency spin-singlet even-parity (ESE), 
even-frequency spin-triplet odd-parity (ETO), odd-frequency spin- 
singlet odd-parity (OSO), and odd-frequency spin-triplet even-parity 
(OTE). 



be performed to search for an odd-frequency superconducting 
state. Though we explicitly studied the N/S junctions only, the 
odd-frequency pairing state is also expected near impurities 
and within Abrikosov vortex cores in even-frequency super- 
conductors. This implies that the odd-frequency pairing is not 
at all a rare situation as was previously considered but should 
be a key concept for understanding the physics of non-uniform 
superconducting systems. 

One of the authors Y.T. expresses his sincerest gratitude to 
discussions with K. Miyake, N. Nagaosa and K. Nagai. This 
work is supported by Grant-in-Aid for Scientific Research 
(Grant No. 17071007, 17071005 and 17340106) from the 
Ministry of Education, Culture, Sports, Science and Technol- 
ogy of Japan. 



and (4)-(7) are worthy of remark. As seen from the above 
Table, the pairing symmetries on both sides of the interface 
are the same, ESE, OSO, ETO and OTE, respectively. As a 
result, the Josephson current can flow across the interface in 
these cases. 

These results can be applied to actual materials. Recently, 
Fuseya et. al. predicted that the OSO state could be realized 
in CeCu2Si2 and CeRhlns [6]. It is consistent with some ex- 
periments 112 ill . Here, we propose a robust check of pairing 
symmetry using the Josephson effect between ESE (conven- 
tional low T c ) and OSO (e.g. CeCoIn 5 ll22ll ) superconductors. 
If two ESE superconductors are attached to opposite (parallel) 
sides of an OSO sample, the ESE order parameters induced 
at the two interfaces in OSO will have opposite signs. Then 
the structure will behave as a 7r-junction. Detection of a tt- 
shift would thus be an unambiguous signature of OSO pairing 
symmetry, similar to the phase-sensitive tests of d-wave sym- 
metry in high T c cuprates [23]. As regards the OTE state, the 
promising system is a diffusive ferromagnet /spin-singlet s- 
wave (DF/S) hybrid structure, where OTE state is induced in 
DF. Recent calculation of the Josephson effect in spin-triplet 
p-wave / DF/S junctions [20] is consistent with the present 
prediction. From this point of view, it is of interest to study 
junctions between Sr9RuOA ll24ll and DF/S hybrids. 

In summary, using the quasiclassical Green's function for- 
malism, we have shown that the odd-frequency pairing state 
is generated near normal metal / even-frequency supercon- 
ductor (N/S) interfaces in the absence of spin flip scattering. 
When the pair potential in the bulk has an even-frequency 
symmetry (spin-singlet even-parity ESE or spin-triplet odd- 
parity ETO state), the resulting order parameter at the inter- 
face has an odd-frequency symmetry (spin-singlet odd-parity 
OSO or spin-triplet even-parity OTE state), in agreement with 
the Pauli principle. On the other hand, if a superconduc- 
tor has an odd-frequency (OSO or OTE) order parameter in 
the bulk, then, respectively, ESE or ETO pairing state should 
be induced near the interface. It follows from the above re- 
sults that the Josephson coupling may occur between odd- and 
even-frequency superconductors and phase-sensitive tests can 
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